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1 Introduction 



In 2 dimensions certain integrable quantum field theories can be restricted to the x < half 
line by imposing suitable boundary conditions without destroying their integrability (1J. In 
addition to their theoretical interest, these models have important physical applications in 
various impurity problems (for a review see @). 

An example is provided by the sine-Gordon model: as was argued in flT], its boundary 
version: 

oo oo 

S= J dt J dxCsc- J dtV B (^ B ), £ 5G = ^(^$) 2 -^(l-cos(/3$)) ) (1.1) 

— oo —oo — oo 

(where $(x,t) is a scalar field, (3 is a real dimensionless coupling and $#(£) = $(x,t)\ x= o) 
preserves the integrability of the bulk if the boundary potential is chosen as 

V b ($b) = -M cos (~($b - <t>o, 

where M and </>o are free parameters. A novel feature of the boundary sine-Gordon model 
(BSG) is the appearance of a complicated spectrum of boundary bound states (BBS) in addi- 
tion to the well known bulk ones HH-fl], f|. The complete spectrum of these bound states and 
a full explanation of all the poles are known only for BSG with Dirichlet boundary conditions 
(which corresponds to taking M — » oo, $e(t) = <p ) Q. 

In this paper we investigate the bound state spectrum in sine-Gordon model with Neumann 
boundary condition (SGN). This boundary condition is in a sense the opposite limit to Dirichlet 
(Mo = 0, thus 0o becoming irrelevant) and has interesting properties as the non conservation 
of topological charge. Furthermore any BSG model, which is not Dirichlet, behaves in the 
ultraviolet limit as if it had Neumann boundary condition. 

We determine the spectrum of boundary bound states by using the bootstrap principle and 
give in fact an inductive self consistent proof of their existence. We also give the explanation 
of all the poles in the various reflection factors in terms of these bound states and Coleman- 
Thun diagrams ||. The importance of Coleman-Thun diagrams in the context of boundary 
bootstrap was first emphasized in 0. An interesting feature of SGN is that we find instances 
when Coleman-Thun diagrams and bound state creation coexist (a similar phenomenon was 
previously observed in the case of boundary Yang-Lee theory [^]). To discover this one has to 
compute and compare the residues of Coleman-Thun diagrams and the reflection factors, thus 
one has to go beyond the usual argument that checks only the existence of a diagram with 
the required pole. 

Once we know the spectrum of SGN on the infinite half line the next step is to inquire 
how it is changing when the model is restricted to a finite line segment < x < L with 
suitable boundary conditions. We develop a transfer matrix formalism and use it to discuss 
the description of boundary bound states in this setting in some detail. Since the resulting 
finite volume spectra depend on the reflection factors, they provide an ideal laboratory to 
check these reflection factors, if we can measure the spectra independently. 

To this end we develop the truncated conformal space approach (TCSA) || to SGN. 
Determining the finite volume spectra numerically by using TCSA we are able to compare 
them to the theoretical predictions, and the complete agreement we find gives a strong evidence 
for the correctness of the various bound states and reflection factors. 
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The paper is organized as follows: in the next section we describe some classical solutions 
of SGN and show that the semi-classical quantization of the boundary breathers explains 
a part of the poles in the breathers reflection factors on non excited boundary. In section 
3, using the bootstrap principle, we derive the spectrum of boundary bound states and the 
associated reflection factors. Section 4 is devoted to the discussion of finite size effects in a 
large volume L, with some non trivial boundary conditions at the ends of L. In section 5, as 
the main step of TCSA, the Hamiltonian of BSG on < x < L is described as that of a bulk 
and boundary perturbed free boson with appropriate boundary conditions. We compare in 
detail the TCSA data and the theoretical predictions in section 6. We make our conclusions 
in section 7. The paper is closed by three appendices: in Appendix A we give in details the 
proof of the existence of the boundary bound states. In Appendix B we explain the poles of 
the soliton and breather reflection factors. Finally, in Appendix C we review those aspects of 
boundary c = 1 theories that are necessary to set up our TCSA. 

2 A few classical solutions and the semiclassical spectrum 
of boundary breathers in SGN 

In boundary sine-Gordon theory with Neumann boundary condition (SGN) the boundary 
potential is absent: M = 0. As a result, <&(x, t) satisfies free boundary condition at x — 
0: d x &(x, t)\ x=0 = 0. This boundary condition preserves the $ <-> — $ charge conjugation 
symmetry (C symmetry) of the bulk theory, but violates the conservation of topological charge. 
This can be seen already on the simplest classical solution of SGN, describing the scattering 
of a classical soliton on the boundary. This solution, $ s (x, t), is obtained by restricting to the 
x < half line the §sa(%, t) soliton anti soliton solution of the bulk theory: 
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sinh (umt / yl — v?) 



$ s (x,t) = &sa( x , t) = — arctan ^=^=- , for — oo < x < 0. (2.1) 

p \-u cosh(mx/v 1 — w 2 )-l 

(Since d x <&sA(x,t)\ x =o = it indeed satisfies the Neumann boundary condition). This solution 
exhibits that the incident soliton reflects as an anti soliton from the Neumann boundary, i.e. 
topological charge changes by two units in this scattering. 

Another observation that plays an important role in the sequel is that the 'standing' 
breather solution of the bulk SG, oscillating with period r around x = 0: 



sin(mt/f ) 



$ r (x,t) = -arctan y/(f ) 2 - 1 — — v / / \^ = , r = — , (2.2) 



(3 L cosh (mx a/ 1 — 



T 



mr 



2tt 



also satisfies d x § T (x, t)\ x= o = 0. Therefore the classical breather solution, bound to the wall 
at x = is given by 

<f> T (x, t) = $ r (x,t), for - oo < x < 0. (2.3) 

In SGN 'boundary dependent' simple poles are found in the reflection factors of the various 
breathers H that may describe boundary bound states (BBS). In the reflection factor of the 
n-th breather, B n , this pole is located at the 6 = inpn/2 [] value of the purely imaginary 

lr The parameter p is determined by the sine-Gordon coupling constant as p = /3 2 /(8n — f3 2 ). 
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rapidity, so that if it really corresponds to a BBS, then the energy of this state above the 
ground state is 

e n — eo = M sm(npn) . (2.4) 

In the form of two lemmas sufficient conditions were given in @ that guarantee that a simple 
pole in a reflection factor cannot be explained by the Coleman-Thun mechanism ||, i.e. that it 
describes a bound state. However, these conditions are not satisfied for any of the poles above, 
thus the question whether they correspond to BBS remains open. Furthermore, a Coleman- 
Thun explanation of a subset of the poles is given by the solitonic version of diagram (c) on 
Fig.( |3.2j ), when the index n of the breather satisfies n > 1/ (2p). (This is discussed in detail in 
the next section). We clarify the status of the poles for n < l/(2p) below by showing that the 
energies associated to them as hypothetical BBS match exactly with the spectrum of bound 
states obtained by semi-classical quantization from the classical boundary breather solutions 
of SGN (|2.2| - |2~lf ). We achieve this by adapting to this problem the semi-classical derivation of 



the breather spectrum in the bulk sine-Gordon theory as given in the classic paper ||. 

The semi-classical (WKB) quantization of any periodic classical solution, (p c i in a field 
theory can be summarized by the equations: 



d_ 



Sd(,<p d ) + Satta) - J2( n i + = E , ( 2 - 5 ) 

i=0 

W {ni} (E) = 2Knh, K integer, (2.6) 



where 



W {ni} {E) = S d {<j>d) + S a te d ) + Er{<j)d) - + (2.7) 

i=0 

Here r(4>d) denotes the period of the classical solution, Vi((f>d) stand for its stability frequencies 
(that characterize the quasi periodicity d(x, t + r) = e tUi £i(x,t) of the solutions of — J^ + J^- — 

(t0O</> 1 £.i( x 't) = 0), and are a set of non negative integers. Furthermore the subscript 
'ct' stands for counter term contributions, whose effect is to cancel divergences in the infinite 
sums. Eq. (|2.5| ) picks out the appropriate classical periodic solution for a given E and a given 



set of integers rij. Then those values of E which satisfy Eq.( |2.6|42T7|) , give the renormalized 
bound state energies. 

With the Neumann boundary condition the boundary potential is absent, and as £sg($t) 
is symmetric for x 1— > —x, we easily obtain 

S d (^r) = ^ farccos (I) - vW 3 ^) , (2-8) 



(3 2 \ \f 

which is just half of the corresponding expression in the bulk theory. The stability frequencies 
for $ T are the same as for $ T , but the fluctuations must also satisfy <9 x £(x, t)\ x=0 = 0. 

Therefore only half of the bulk stability fluctuations appear in the Neumann problem, namely 
those that are even under x 1— > —x, (and the same applies to the fluctuations contributing 
to the vacuum energy E vac appearing in the counter term). Thus effectively we must take 
Yl &i = \ v%. To obtain the basic boundary bound states we set all rii = 0, and using the 
explicit form of the counter-terms and the sums over Vi as given in [|] we find finally 

S a ($ T )-J2h = -^Sci(®T). (2.9) 
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Using these results in Eq. (|2.5| , |2.7|) leads to 



E = — , W(E) = -arcsin — — . 2.10 

np r p \ m J 

As a result the quantization condition, W[Ek) = 2ttK, gives 

E K = Msm(KpTr), (2.11) 

where M = ^ is the semi-classical soliton mass. Clearly this reproduces (|2.4|) . Since t(Ek) ~ 
^(e'ptt) the range of if-s where bound states exist (i.e where t(E k ) < oo) is K < l/(2p). 



3 Boundary bound state spectrum from bootstrap princi- 
ple 

In this Section we determine the spectrum of boundary excitations, the related soliton and 
breather reflection factors and show how their poles can be explained in terms of on-shell dia- 
grams. We start with a summary of the bulk scattering properties and then review the result 
of Ghoshal and Zamolodchikov concerning the reflection factors of the soliton and breathers on 
the non-excited boundary. From their pole structure and the bootstrap equation we conjecture 
the minimal spectrum of the excited boundary states (or with other words boundary bound 
states), and prove in Appendix A that they can be created by successive soliton absorptions 
on the wall at purely imaginary rapidity. Having determined the reflection factors we explain 
all their poles in terms of on-shell diagrams which correspond either to some Coleman-Thun 
mechanism or creation of boundary bound states. In this study we rely heavily on the ma- 
chinery worked out in |4] for the Dirichlet case. However, in contrast to the Dirichlet case we 
find instances when a Coleman-Thun type diagram happens to give a pole of the same order 
as the reflection amplitude has, but with wrong residue, thus leaving the possibility of excited 
boundary state creation at the same time. 



3.1 Bulk scattering properties 

In the bulk SG model any scattering amplitude among solitons and anti solitons factorizes 



into a product of two particle scattering amplitudes, of which the independent ones are [|T0 

r(2(/- i)A-^)r(2/A + i 
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Since we are concentrating on the bound state poles located at purely imaginary rapidities we 
use the variable u instead of 6 and refer to it as the rapidity from now on in this Section. The 
other scattering amplitudes can be described in terms of the functions 



{y} 



m (#) 
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(a) Breather fusion 



(b) Soliton and anti soliton 
fuse to a breather 



Figure 3.1: Bulk vertices 



as follows. For the scattering of the breathers B n and B m with n > m and relative rapidity u 



we have 10 



S nm {u) = S™{u) = {n + m- l}{n + m - 3} . . . {n - m + 3}{n - m + 1} , 
while for the scattering of the soliton (anti soliton) and B n we have 

{1 + A} if n is even 



fi» = S+M = ST>) = {n - 1 + \}{n - 3 + A} . . . 



\/{A} if n is odd 



All the poles of the scattering amplitudes in the physical strip originate from virtual processes 
either in the forward or in the cross channel of the diagrams (a-b) on Fig. (|3.1|) , where the 
useful definition 



u,. 



nir 
2A 



was also introduced.^ For each such process a coupling as /™^ m or can be attributed and 
it is known that = (-l) n /™ + . 



3.2 Reflection factors on ground state boundary 

The reflection factor of the soliton on the ground state Neumann boundary was found by 
Ghoshal and Zamolodchikov 0. For the topological charge preserving process, (i.e. when a 
soliton (anti soliton) reflects as soliton (anti soliton)), we have 



P(u) = cos(\u)Ro(u)a ^(A + 1), uj tr(0, u) 



where 



Ro(u) 



n 

i=i 



r(4/A-^)r(4A(/-i) + i-^) ;/ 

r((4/-3)A-^)r((4/-l)A + l-^) A ; 



2 On all the space-time diagrams time develops from top to bottom. Solitons (or anti solitons) are denoted 
by solid lines, while breathers by dashed ones. 
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is the boundary condition independent part and 

r(i + a + (21 - i)A - £ )r(i - % + (21 - i)A - ^; 



, N COS X 

a(x, u 



cos(x + \u) 



n 

z=l 



r(| - a + (2/ - 2)A - &)T{\ + £ + 2ZA - 



describes the boundary condition dependence.(p?he poles of a(x,u) in the physical strip are 
located at u = | — «2fc+i or at 7r — | + W2fc+i, > ; and it has no zero there. 

For the topological charge changing process, (when the soliton (anti soliton) comes back 
as anti soliton (soliton)), we have 

sin(Au) 
sm(^j-) 

Note that the topological charge is changed by two in this process, thus the parity of the 
soliton number is conserved and we have an odd and an even sector. 

The breather reflection factors share the same structure as the solitonic ones as a conse- 
quence of the general expression in ||. For Neumann boundary condition and breather B n 
they have the following form 



R^(u) = Rf\u)S^ (|(A + S^(0,u) 



where 



nWf.A - (2) (2A + X ) TT ttx) (L + X ) . q(n) (r \ _ TT (xtt 2 + a) 

\2X ~ 2) l=\ \2X ' 2/ ;_l-n \ Xtt ~ 2 ~ X) 

In general would describe the boundary independent properties and the other factors 
would give the boundary dependent ones. In the Neumann case, however, there are coinci- 
dences among the poles and zeros of the various factors and it is better to rewrite the reflection 
factor as 



R (n) (u) = (2n-l)(2n-3 + 4A)...(2(n-A;)-l + (l-(-l) fe )2A)...| ^ 
where 

(y±±\ (y-i\ 

(y) ~ IA IA 



(1) if n is odd 
■f 4A) if n is even 



(y+1 _ i\ ( y- 1 _l n 

I 4A 2M 4A T 2j 

is the "half-block" with respect to {y}. 
3.2.1 Pole analysis 

Now consider the poles of the soliton reflection factor P(u). 



• There are boundary independent poles in the physical strip due to the factor Ro(u). 
They are located at u n , n = 1, 2, . . . from which every odd is canceled by the zero of 
cos(Aii). The remaining poles can be described by diagram (a) on Fig.(|3.2|). This is 
consistent with the fact that the C symmetric Neumann boundary can emit and absorb 
only even breathers with zero momentum. 

3 The function a(x, u) has among others the property a(x, u)a(x, —u) — cos ( x+ xu) cos{x-Xu) which corrects 
a typo in and §J. 
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(a) Soliton bulk pole (b) Breather (c) Breather emis- (d) 

triangle son reflection bound 

state 
replace- 
ment 



Figure 3.2: A few boundary diagrams 



The boundary dependent poles in the physical strip are located at 



7T 



U 2r 



71 = 0,1,..- 



The pole at | is responsible for the emission and absorption of a zero momentum soliton 
or anti soliton. For each of the other poles we associate a boundary bound state with 
energy 

e n - e = M sin u 2n 

From now on we put the boundary ground state energy to zero eo = for simplicity. 



Clearly the basically bulk process (a) on Fig. (|3.2j) does not exist for the reflection amplitude 
Q(u) and the factor sin(Aw) takes care of this. The boundary dependent poles of Q{u) coincide 
with and can be explained in the same way as those of P(u). 

Now let us focus on the poles of the breather reflection factors. As a warm-up example 
consider the reflection factor of the third breather: 



R( 3 )(u) = (5)(3 + 4A)(l) 



• It has a simple bulk pole at | — w 3 whose explanation depends on the value of the 
parameter A. If B 6 is in the spectrum then we can use the 'all-breather' version of 
diagram (a) . On this we mean by replacing each soliton line by a B s line having rapidity 
| — u% and B 6 is formed at the vertex. Now if B 6 is not in the spectrum then the 
soliton version of diagram (b) applies. In this case the breather triangle is turned into a 
soliton-anti soliton triangle. 

• The pole at ~ — U\ is of second order and can be explained by diagram (c). 

• The pole at u 3 can be explained by the soliton version of diagram (c) for 3 > [~] 
(which is necessary to guarantee that the soliton in the middle also travels towards the 
wall). This diagram gives a second order pole and the reflection factor has no zero 
at —V3 — uq — f . We have an analogous diagram, however, by changing the amplitude 
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P(u) for Q{u). They differ by a factor of sin(— Az/ 3 )/ sin(4^) = — (— l) 3 but we know that 
/+_ = (— 1) 3 /-+- This gives the sign difference which is responsible for the Coleman- 
Thun type cancellation in the sum of the two diagrams. In the 3 < [§] case we have a 
boundary bound state with energy 

e 3 = M3COSM3 = 2M sin u 3 cos u 3 = Msinw 6 , 

which is exactly the same as the boundary bound state created by the soliton. The 
identity of these two boundary bound states can also be confirmed by showing the 
agreement of the other higher spin conserved charges. 

• For the other two poles at U\ and u 2 we have diagram (b). In the case of u 2 the second 
breather hits the boundary while in the case of Ui the first. 

We could be satisfied now and turn to the general case. There is a subtlety, however, since we 
have not computed the residues. Performing the computation reveals that the residues of the 
reflection factors and diagrams agree except for the pole u\ where they have different signs. 
Thus we must assume that a boundary bound state exists with energy M 3 coswi = ei + e 2 . We 
note that in the Neumann case the bulk pole of Rf\u) and a boundary pole of (|(A+1), u) 
at Mi coincide but a zero in the boundary factor reduces the singularity to first order. If 
we move a bit away from the Neumann point, when the boundary dependent part becomes 
S^\t;(\ + 1) — 5\,u), the two poles split and the bulk pole at u\ can be described by the 
diagram above, while the pole at U\ — 5 corresponds to a boundary bound state. Computing 
the residue and taking the 5 — > limit do not commute, however, resulting in a sign difference. 
This sign difference persists when the state with energy e 2 is not in the spectrum i.e. for 114 > |. 
In this case we have diagram (d), where B 3 is decaying into a soliton and an anti soliton. Then 
the soliton creates the state with energy t\ on which the other anti soliton is reflecting with 
rapidity 114 — |. This diagram is naively second order, however, the same Coleman-Thun 
mechanism that we used in the previous example, takes care of the cancellation. 

From this example we can draw two conclusions. First, breathers and solitons can create 
the same states. Second, even if a diagram describes a pole of the same order as the reflection 
factor has, the residues do not always coincide, thus allowing the creation of a boundary bound 
state in the same time. 

Now we are ready to discuss the general case. The reflection factor R^ n \u) has the following 
poles in the physical strip. 

• There is a simple pole at | — u n . It can be explained in terms of the breather version of 
diagram (a) by forming B 2n or if it is not in the spectrum then by the soliton version of 
diagram (b). 

• The double poles at | — u n - 2 k can be explained by diagram (c). 

• For even breathers there is an extra simple pole at |. It corresponds to the emission of 
a zero momentum even breather. 

• The simple pole at u n corresponds to the creation of the boundary bound state with 
energy e n = Msinw 2n , if it is in the spectrum. Alternatively, it can be explained in 
terms of the solitonic version of diagram (c) where, just as in the example, a Coleman- 
Thun type mechanism ensures the pole to be simple. 
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• If n + k is even then the pole Uk is responsible for the creation of a boundary bound state 
with energy e^+k + e»-t , if e^+k is in the spectrum. Alternatively, we have diagram (d) 
where the state with energy is created in the middle and the Coleman-Thun argument 
has to be used for the scattering. 

• If, however, n + k is odd then the pole at Uk has its explanation in terms of diagram (b). 



3.3 Excited state reflection factors 

We have seen that both the pole of the soliton reflection factor at v n and the pole of the 
reflection factor of B n at u n give rise to a state with energy e n . The identity of these two states 
can be further confirmed by showing that any reflection factors on them agree. Computing 
these reflection factors from the solitonic and breather bootstrap equations gives identical 
results. Therefore we denote these common states by |n). 

In computing the soliton excited state reflection factor we use the 

P\ n) {u) = P(u)S n (u + u n )S n (u - u n ) 

"soliton-breather" bootstrap equation. The result is 

P\n)(u) = a n (u)P{u) ; a n (u) = {2n - 1 + \}{2n - 3 + A} ... {1 + A} . 

We have checked that the purely solitonic bootstrap equation gives the same result. Since 
the solitonic reflection factor has one extra pole at z/ = | we can bootstrap on this pole and 
obtain P\o)(u) = P{u). This means that the ground state is non degenerate, so the creation 
or annihilation of a soliton or anti soliton does not change the boundary ground state. (We 
have already used this property in the explanation of the ground state breather poles) . The 
existence of a unique ground state implies that C symmetry is not broken. Note, however, 
that this is strikingly different from the case for Dirichlet boundary condition at £ = |(A + 1), 



where one has two vacua and the Z 2 symmetry is spontaneously broken pi 



Repeating the same bootstrap program for the amplitude Q{u) we arrive at 

Q\ n )(u) = a n (u)Q(u) . 

For the breather reflection factors on the excited state \n) the bootstrap equation takes the 
form 

R§\u) = R (m \u)S nm (u + u n )S r 
From here we have explicitly 



it - u r 



|n) 



U 



where 



. {to + 2n - 1 - 21} . . . 



b™(u) = {m + 2n - l}{m + 2n - 3} . . . 

{m - 2n + 1} if to - 2n + 1 > 
{2n-TO+l} if 2n - to + 1 > 

Let us focus on the solitonic reflection factors P\ n )(u). The prefactor a n {u) has simple poles 
at vq and v n and double poles at Vk > k = 1, 2, . . . , n — 1. The locations of these poles are 
the same as that of P(u), thus only the order has increased for k < n. We conjecture, and 
in fact prove in Appendix A, that the other simple poles at Uk , k > n, which do not come 
from a n (u), are responsible for the creation of the excited boundary bound states with energy 
e n + efc which we denote by \n, k). 
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Figure 3.3: Various ways of creating the state \n, k). 



Now it is interesting to use the integrability of the model. Using the factorizability argu- 
ment we can shift the soliton trajectory and obtain the three diagrams on Fig. (|3.3|) . 

In the first picture the soliton creates the boundary bound states |n) at an angle v n then 
an other soliton creates \n, k) at an angle v k where k > n. They scatter on each other at an 
angle U2(k-n) which is not proper for forming a breather. Now if we move the second soliton 
trajectory up we obtain the picture in the middle. Here the soliton at v k creates the state \k) 
then the other soliton reflects on it at an angle v n . This reflection factor has a third order 
pole for which diagram (a) on FiR. flB.3l ) shown in Appendix B can be associated showing 
the creation of the state \n, k) as it was expected from consistency. In the third picture the 
excited boundary \k) decays by emitting a soliton. This soliton then scatters on the other 
soliton forming B n+k with rapidity u k - n which finally creates the state \n, k). 

Now we can use any of the processes above (or their breather versions) to compute the 
reflection amplitudes. They turn out to be 

P\n,k)(u) = a n {u)a k {u)P{u) ; Q\ n ,k)(u) = a n (u)a k (u)Q(u) 

and 

From this it is clear how it goes on, so we turn to the description of the spectrum of the 
boundary bound states and their reflection factors in general. 

3.4 The spectrum of the boundary bound states and the associated 
reflection factors 

From the previous discussion it is clear that a boundary bound state can be labelled by a 
sequence of increasing positive integers and denoted by |n l5 n 2 , . . . , n k ) where > n^and 
n k < [|] . Its energy can be written as 

k k 

m\ nun2) ... )nk ) = smu 2n] = y^e n? .. 

The first task is to prove the existence of these states. Unfortunately the proof of |4|], which 
uses the breathers to create the boundary bound states, does not apply. This proof is based 
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on a special property of the breather reflection factors with Dirichlet b.c.-s, namely that to 
every pole either a bound state or a Coleman-Thun diagram can be associated. In contrast 
to the Dirichlet case we have shown that Coleman-Thun diagrams and boundary bound state 
creations may coexist in the Neumann problem. We prove, however, in Appendix A that soli- 
tons at reflection angle v m on the state \rii, n 2 , . . . , n k ) can create the state \ni,n 2 , ■ ■ ■ , n k , m) 
if m > rik- Now the next step is to compute the reflection factors and explain all the poles in 
terms of diagrams. The computation is quite straightforward and the result is 

P\m,n 2 ,...,n k ) (u) = a ni {u)a n2 (u) . . . a nk (u)P(u) 
Q\ ni ,n 2 ,...,n k )(u) = a ni (u)a n2 (u) . . . a nk (u)Q(u) 

and 

<L..,n fc) W = TOO) • ..bZ{u)R {m \u). 

Using the factorizability argument leads to the following creation-annihilation rules: For the 
soliton with rapidity v m 

|ni, . . . , Hi, n i+1 , ...,n k ) — > \ni, . . . , n*, to, n i+1 , ...,n k ) 

if rii < m < n i+ i. For B m with rapidity U2 ni + m we have 

|ni, . . . , Hi, . . . , n k ) — > |ni, . . . , rij + m, . . . , n k ) , 

while for B m+l with rapidity u\_ m and / > m the process is 

|ni, . . . , rii, n i+ i, . . . , rij, n j+1 , ...,n k ) — > \ni, ...,rii,m, n i+1 , ...,nj,l, n j+1 , ...,n k ) ■ 

Now the most difficult problem is to explain all the poles of the reflection factors in terms of 
diagrams. Since it is quite cumbersome we put it in Appendix B. 



4 Finite size effects in large volume 

4.1 Bethe-Yang equations for particles in finite volume with bound- 
aries 

Imagine putting N particles in finite volume L, such that the ith particle Ai has mass Mj 
and rapidity (which we choose nonnegative since only the absolute value matters, due to 
reflection on the boundaries), plus some additional internal quantum numbers ctj (distinguish- 
ing it within a multiplet of the same mass Mj). Let us denote the reflection factor of the 
ith particle on the left/right end of the interval by R^(9i) and R^(9i), respectively and the 
scattering matrix of particles Ai and Aj by S^\9i — 9j). Obviously in the general case these 
are matrices with the following structure 

R R^i) = j-RflVOSj > 

S^\9 i -9 j ) = {s^\9 i -9 j )Tc?\ ■ 
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We define a family of transfer matrices that acts in V\ <g> . . . <E> Vjv, where Vi is the internal space 
labelled by the multiplicity index «j of the ith particle. The transfer matrix T k describes a 
unitary transformation on the wave function resulting from the scattering of A k on the other 
particles and the two boundaries when we take the particle around the volume (starting e.g 
to the left, reflecting on the left boundary, coming back to the right, reflecting on the right 
boundary and returning to the original position): 

t£'"' n) (e 1 , ...,e N ) = R%\e k ) n s iJ ' k \e k + ej)R^\e k ) J] s^ k \e k - o,) . 

We omitted the multiplicity indices of the matrices (these are indicated by the upper indices 
in brackets and can be restored easily). The scatterings were performed in a particular order, 
due to the equations of factorized bulk and boundary scatterings, the order is eventually 
immaterial: performing the scatterings in any other way results in the same transfer matrix. 
Another consequence of factorized scattering is that the matrices T k form a commuting family 
for a given set of rapidities. The wave function can be characterized by a vector ip ai ...a N living 
in the tensor product space, and transforming under all such monodromy operations as: 

exp(2iM k Lsmh(6 k ))Tj: 1 - N) (6 u ...,6 N )^ = ij , k = l,...,N. (4.1) 

The prefactor in this equation is the phase acquired by A k due to its momentum (whose 
absolute value is M k smh9 k ). Equations (|4.1|) are the so-called Bethe-Yang equations in the 
presence of a boundary. 

Since the transfer matrices commute, they can be diagonalized simultaneously. Let us 
denote the eigenvalues of the transfer matrix T^ 1 ''"'^ . . . , 9 N ) by (9i, ... , 9^), where 
s = 1, . . . , D N enumerates the eigenvalues (with multiplicities) and D N = dim Vi <g> . . . ® V^. 
The corresponding common eigenvectors will be called tp( s > (9i, . . . , 9n)- 

The solutions to the Bethe-Yang equations ( ^.1| ) are given by the wave function amplitudes 

provided that the rapidities solve the equations^ 

exp (2zM fe Lsinh(0 fe )) (9 U ...,6 N ) = l , k = l,...,N (4.2) 

for some (fixed) s. These equations effectively describe the quantization of momentum in finite 
volume with boundaries. The total kinetic energy of the particles is then given by 

N 

Ekin = Mj cosh 9j , 

3=1 

which is the energy difference with respect to the state with no particles in this approximation. 

The above calculation supposes that the particles do not overlap with each other and the 
boundary substantially, i.e. it is valid only when the volume is large enough: MjL ^> 1 for all 
i. Note that if the particles involved in the calculation have no multiplicity labels, the transfer 
matrix is a scalar phase itself and can be directly substituted into Eqs. (|4.2j ). Taking the 
logarithm of Eqs. (|4.2[) yields 

2M k Lsmh(e k )-ilogX k s) (9 1 ,...,9 N ) = 2irl k , fc = l,...,JV. (4.3) 
4 For scalar particles, these equations appeared in Eg]. 
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where the Ik are integers that we call Bethe quantum numbers. In this way a given multi- 
particle state can be labeled by giving s and a set of integers I\ , . . . , In and this labeling 
is normally unique given a consistent choice of the branch of the logarithm. There are some 



exceptions when for given s and I\ , . . . , ijv Eqs. (4.3) admit more than one solution; such 
cases will be explicitly indicated below. 

Note that in these considerations we neglected the contribution of vacuum polarization. 
Since, on dimensional grounds, we expect this contribution to behave as C(exp(— ML)) for 
large L-s (with M being some characteristic mass scale), the polynomially decreasing, leading 
finite size corrections are indeed described by Eq.(|4.3|). 

4.2 Some particular cases 

Here we write down and analyze Eqs. (|4.3| ) for some simple configurations that will play an 
important role in the sequel. 

4.2.1 States containing a single scalar particle 



For a state containing a single scalar particle of mass M, Eqs. ( |4.3|) reduce to the single 
equation 

2ML sinh(#) - i log R L (6) - i log R R (6) = 2nl . (4.4) 
The energy of the state with respect to the vacuum can be written as 

E(L) - E (L) = E R + E l + Mcoshfl , 

where E Ri l denote the boundary energy of the left and right boundaries with respect to the 
boundary ground state (note that the boundary can also be in an excited state). One can 
easily express the function E(L) — E (L) in a parametric form using 6. 

Our studies show that for ML <C 1 Eqn. (|4.4|) always has a real solution irrespective of 
the choice of the quantum number /. However this may not be true for larger values of ML. 
Let us discuss the particular case 1 = and introduce the notation 

d 

Pl,r= i-^logR LtR (9) 

It is a generic property of reflection factors that zlog Rl,r(0) has a finite limit as 9 — > ±oo; 
it also follow from unitarity that they are odd functions. Therefore the equation 

2 ML sinh(#) = i log R L (6) + i log R R (6) (4.5) 

has a pair of nonzero solutions as long as 

2ML <p L + p R . 

Increasing L, these two solutions move towards the origin, which they reach at a finite volume 

Lo = ^(Pl + Pr) ■ 

The question arises: what happens for volumes L > L ? If ML ^> 1 the Bethe- Yang 
equations should give a good approximation of the finite volume spectrum around this point. 
However, eigenvalues cannot simply disappear as the Hamiltonian we are concerned with is 
Hermitian. In fact, the spectrum is real and every eigenvalue is a continuous (even smooth) 
function of L. This leads us to the description of boundary bound states. 
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4.2.2 Description of boundary bound states in finite volume 

Quite naturally, the two solutions of Eqn. (|4.5|) do not disappear, but continue their life as 
complex (eventually imaginary) solutions. Substituting 6 = iu into Eqn. ( |4.5|) we obtain 

2ML sin(w) = log R L (iu) + log R R (iu) . (4.6) 

As reflection factors always take real values on the imaginary axis (as a result of analytic 
S'-matrix theory) this is a real equation and it has a pair of real solutions for u exactly when 
L > L . The energy corresponding to this solution is 

E(L) - E (L) = E R + E L + Mcosu . 

This is smaller than E R + E L + M which means that this state cannot correspond to a real 
particle between the two boundaries. 

The reflection factors Rl and R R normally have poles on the imaginary axis. Let us denote 
the one which is closest to the origin by u* and suppose it occurs in Rl onlyQ. Then one has 

u — > u* as L — > oo 

and so 

E(L = oo) - E (L = oo) = E R + E l + Mcosu* . 

Suppose also that u* corresponds to a boundary bound state \B* L ) in the bootstrap, then the 
energy of such a state is 

E L * = E L + Mcosu* . 

It is then tempting to interpret the corresponding state as one without particles but in which 
the left boundary is excited to the state corresponding to the pole at « = u* . Indeed our 
numerical data show complete consistency with this interpretation. 

When the left and the right boundary conditions are identical (e.g both are Neumann), 
both reflection factors have a pole at u = u*. One can then interpret the resulting 'finite size 
state' as one in which one of the boundaries is an excited one. This can be realized with two 
wave functions (as we have two identical boundaries) and we expect that the state described 
above corresponds to one of them. Note that in the above solution we always have u < u*. 
The other one can be described (at least for large enough L) by noting that there is going 
to be another solution u' to Eqn. ( [4.6| ) which satisfies u' > u* and also approaches u* for 
L — > oo. It is clear, that, for finite L, the solution with u' has a lower energy than that of 
with u. 

When L and R are not identical and the boundary state, \B R ), is not identical to the 
excited \B* L ) either, one again expects two wave functions for the system with no particles, 
so the two solutions that exist in this case too can be interpreted once again in this way. 
However, when \B R ) = \B* L ) only one of the solutions can have an interpretation in terms of a 
real physical state. This shows that not necessarily all solutions of the Bethe-Yang equations 
(|4~1~1) correspond to physical states. In fact, when u* corresponds to a Coleman-Thun pole, we 
do not even expect the 1 = state to appear in the spectrum, since then we would have some 
state for L < L that does not have any physical continuation to L < L. 

There is a generalization of this line of thought to other poles that are further away from 
the origin in the physical strip < u < n/2 . We cannot expect these to be realized in 

5 We meet this situation e.g. in sine-Gordon model with <f)\n — (special Dirichlet) and d x (f>\L = 
(Neumann) b.c-s. 



15 



an analytic continuation from a state with real rapidities, since any such continuation goes 
through the origin u = 9 = and necessarily runs into the pole closest to 0. However, we 
can calculate solutions of Eqn. ( |4.6| ) around any such pole (of course one cannot go too far 
from the given pole since one then generally runs into the domain of some other singularity of 
the reflection factor). For large enough L, such a solution exists and generalizing the above 
scheme one expects that it can describe the leading finite size behaviour of some excited 
boundary state in finite volume if the pole in question participates in the bootstrap. Once 
again, Coleman-Thun poles can have no such corresponding state in the finite size spectrum. 

A given excited boundary state may arise from different reflection factors as a pole (e.g. in 
SGN the state |n) arises from both the soliton and the B n reflections) and therefore a further 
possibility is that it can be described in finite volume in more than one ways. 

Later in this paper we show that all the above considerations are well confirmed by our 
numerical studies. We wish to remark that similar ideas to describe bound states appeared in 



13 for the bulk and in 14] for the boundary case 



4.2.3 States with two scalar particles 

In the case of two scalar particles of mass Mi and M 2 , the Bethe-Yang equations take the 
form 

2M X L sinh 9 1 - i log log R$ (0i) - i log S {12) (9 1 - 9 2 ) 

2M 2 L sinh 6 2 - i log R {2) (6 2 )-i log R® (0 2 ) - % log S {12) (6 2 - t ) ' 

-zlog^ 12 )(0 2 + ^i) =2tt/ 2 

E(L) - E {L) = E R + E L + M l cosh0! + M 2 cosh0 2 • 

In this case, these equations can only be solved by numerical iteration. One can also continue 
these equations to imaginary value of one (or both) of the rapidities. Using that the reflection 
factors are real for purely imaginary rapidities together with unitarity and real analyticity of 
the bulk S matrix one can show (just as for the one-particle case) that the resulting equations 
are consistent and the energy of the solutions is always real. Consider e.g. the Bethe Yang 
equations for B n -B m with rapidities iu and 6 2 respectively, which, before taking the logarithm, 
can be written as 

R < £ ) (iu)R < S ) (iu)\S M (9 2 + iu)\ 2 e~ 2m " lsinu =1 , 1 = ML, 
DM//, \pM/ fl \ S {n ' m) (6 2 + iu) i2mmlsinhd2 _ M ra , m (4-8) 

The first is a real equation with real entries, and each factor on the left hand side of the second 
equation is of unit modulus. Thus, when taking the logarithms, for the complete system there 
is only one quantum number, J m , coming from the second equation. Two particle states with 
one rapidity being imaginary are important, as one can argue that the best way to describe 
one particle states moving between an excited and a ground state boundaries is to use ( |4.7]44.^) 
with ground state reflection factors and an imaginary rapidity tuned into the vicinity of a pole 
corresponding to the excited state. When both rapidities are imaginary, the solutions of ( |4.7|) 
describe zero particle states with excited boundaries. 
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5 Truncated Conformal Space Approach (TCSA) for the 
boundary sine-Gordon model 

Here we describe the Hamiltonian of boundary sine-Gordon model (BSG) living on the line 
segment < x < L as that of a bulk and boundary perturbed free boson with suitable 
boundary conditions. This is the starting point of the TCSA analysis. 

The basic idea of TCSA is to describe certain 2d models as relevant perturbations of their 
ultraviolet limiting CFT-s ||. If we consider boundary field theories, then the CFT-s in the 
ultraviolet are in fact boundary CFT-s. The use of TCSA to investigate boundary theories 
was advocated in [|14], [T5|| . 

As the bulk SG can be successfully described as a perturbation of the c = 1 free boson 



lGj, it is natural to expect, that the various BSG models are appropriate perturbations of 
c = 1 theories with Neumann or Dirichlet boundary conditions. Therefore we take the strip 
< x < L and consider the following perturbations of the models described in detail in 
Appendix C: 

poo i>L 



S = J (^dp$d»$ + iicoa(J3&)^dxdt + 



"OO 

+ / (M, 



COS ( ^($ B - <f> ) ) - 1 



+ M L 



dt . 



Here, for finite M's, Neumann boundary conditions are imposed in the underlying c = 1 theory 
on the boundaries, while if any of the M -s is infinite then the corresponding term is absent 
and the boundary condition in the underlying conformal theory on that boundary is Dirichlet. 
We can rewrite the Hamiltonian of the system in terms of the variables associated to the plane 
using the map (x, it) = £ — > z = e 1 ^, and by changing the integration variable we have 

H =H CFT + | Q 2/1/5-1 jT (V p {e», e' ie ) + V. p {e*, e"")) d6+ 
ilT + ^*_| (1)) + 

Mr /1T\hp 



^F(l) (^*f(-i) + ^*_ f (-i) 



Now the computation of the matrix elements of the bulk and boundary vertex operators V±p 
and \&±/?/2 (with conformal dimension hp) between the vectors of the appropriate conform 
Hilbert spaces is straightforward. Also the integrals can be calculated explicitly. Truncating 
the Hilbert space at a certain conformal energy level E cnt (which is nothing but the eigenvalue 
of the zeroth Virasoro generator) and diagonalizing the Hamiltonian numerically we arrive at 
the TCSA method. 

The TCSA Hamiltonian for BSG with Neumann boundary conditions at both ends is 
obtained from ( |5.1| ) by setting M = M L = and using ( |C.5| ), ( |C.6| ) with n = 1, r = V^r / ' (3 
for He ft and Vp, while the TCSA Hamiltonian for BSG with mixed boundary conditions is 
obtained by setting M = Ml = and using QC.9D and flC.lOD with the same n and r for 
He ft and Vp. 

We choose our units in terms of the soliton mass M. The bulk coupling fi is related to M 

by 

m = «03)m 2 - 2 ^, fy» = ir. ( 5 - 2 ) 
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where re(/3) is a dimensionless constant. In the bulk SG, from TBA considerations, the exact 
form of k((3) was obtained in ||17|, and we use the same form also here in BSC Once we 
expressed /i, the Hamiltonian, ( |5.1| ), can be made dimensionless h = H/M, depending only 
on the dimensionless volume parameter / = ML and (3 (in the general case it also depends on 
M 0iL /M and O , L ). 

6 TCSA analysis of SGN in finite volume 

In this section we review the TCSA analysis of the BSG with Neumann boundary condition. 
The aim of this investigation is twofold: first we want to verify the various reflection factors 
and the boundary spectrum obtained from the bootstrap procedure, and second we want to 
obtain further information about the finite volume behaviour of SGN. We consider in some 
detail the questions of the ground state and the associated boundary energies, the low lying 
one particle states and their reflections on the ground state wall and finally the new states 
predicted by the bootstrap. 

In SGN, as C symmetry persists, there are two sectors, the C = 1 even and the C = — 1 
odd ones. The bulk breathers naturally belong to one of them, as the C parity of the n-th 
breather is (— l) n . However, since solitons and anti solitons can reflect into themselves as well 
as into their charge conjugate partners, solitonic one particle states (i.e. states, whose energy 
and momentum are related by E = \J P 2 + M 2 where M is the soliton mass) are there in both 
sectors. The ground state is in the even sector, while we expect |1) to be of the lowest energy 
one in the odd sector (at least for p < 1/2, which is true for all the cases investigated by us). 

If at both ends of the strip Neumann boundary condition is imposed (NN case), then 
sometimes we have to face the problem of identifying the TCSA lines (states) corresponding 
to the symmetric/antisymmetric combinations of some identical excited states of the identical 
boundaries, as described in section 4. Therefore in the NN case the finite volume spectrum and 
consequently the TCSA line sequences are more complex than the 'naive' one derived by the 
bootstrap procedure in section 3. To avoid this complication when 'verifying' the Neumann 
spectrum we consider a system, where Dirichlet boundary condition with 0o — is imposed 
on one of the boundaries, while keeping the Neumann one at the other (DN case). Since this 
mixed boundary condition is also C symmetric we do not lose the presence of even and odd 
sectors on the one hand, while, on the other, we can be sure, that any boundary bound state 
(BBS) found in this case can be attributed to the Neumann end of the strip, as this Dirichlet 
boundary condition has no bound states Q. 

As charge conjugation acts on the fundamental scalar field by $ i— ► — $, it is straightforward 
to implement the projection onto the even and odd sectors in the conformal Hilbert spaces 
used in TCSA. This projection has two beneficial effects: on the one hand it effectively halves 
the number of states below E cut , thus it drastically reduces the time needed to obtain the 
complete TCSA spectrum, and on the other the separate spectra of the even and odd sectors 
are less complex and therefore easier to study than the combined one. 

In our numerical studies E cut varied between 16 and 25, and this resulted in 4.5 x 10 3 
- 8 x 10 3 conformal states per sectors. In the DN case the number of states below E cut is 
independent of p, while in the NN case it depends sensitively on it. 
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Figure 6.1: Ground state energy versus I in the DN case at various p-s. The slope of the 
dashed lines is — |tan(^ E ). 

6.1 Ground state and boundary energy 

In any boundary field theory in finite volume, for large enough L-s, we expect on general 
physical grounds, that the volume dependence of the ground state energy can be written as 



where e is the ground state energy density and £l,r are the ground state boundary energies 
associated to the left and right ends of the strip. In the BSG model one can argue that e 
is nothing but the well known bulk energy constant, — ^f-tan(^), of the bulk SG theory. In 
BSG with Dirichlet boundary conditions, LeClair et al. (1^| were able to obtain <?Dir from 
TBA, at least for those values of p when not only the boundary but also the bulk scatterings 
become diagonal. The case of Neumann boundary condition is more complicated as there is 
no p where the boundary scatterings would become diagonal. 

Nevertheless using TCSA, we can get information on ^Neu- To this end we determined 
Eq(1)/M (M is the soliton mass and I = ML is the dimensionless volume parameter) from 
our TCSA data at various values of p both in the (DN) and in the (NN) cases; the results are 
shown on Figs. (|Q| - |6.2| ) These figures show in a convincing way that in these BSG models 
e is indeed the expected one. Interestingly the (DN) data for E (l)/M are consistent with 
a straight line passing through the origin; i.e. with S^ eu (p) = — £vir(p, <fio — 0)- Since from 
ref.[l8|] we know £oir analytically, using the (DN) data we have a prediction for E (l)/M in 
the (NN) case: a straight line with slope parameter — |tan(^) and intersection 



E {L)=e L + £ L + £ R 



(6.1) 




(6.2) 



Fig. (|6.2j) shows that this prediction is consistent with our (NN) data. 
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Figure 6.2: Ground state energy versus I in the NN case at various p-s. The slope of the 
dashed lines is the same as on Fig. (|6.1|) , but their intersection is given by Eq.(|0|). 



6.2 Low lying one particle lines, reflections on the ground state wall 

Next we summarize what we can learn about the Neumann spectrum by studying the lowest 
lying one particle states. 

6.2.1 Breather lines 

Using the formalism developed in section 4 it is straightforward to express in parametric form 
the Bethe - Yang (BY) lines corresponding to the various breathers moving between ground 
state walls at the ends of the stripQ : 

(E(L),L) = (M n cosh6 1 (2TTl n + t\ogR ( £\e) +i\ogR^ ) (6))/(2M n smhe)) . (6.3) 

Here E(L) is the energy of the n-th breather above the ground state, 6 (the parameter of the 
line) is the rapidity of B n (which, therefore, necessarily must be real and non negative), and 
M n is the mass of B n . I n is the quantum number characterizing the line; if R^^R^ (0) = 1 
then I n is integer, while if R^ ] (0)^(0) = -1 then I n is half integer. (If i*[ n) (0) = R { r\0) as 
in the (NN) case, then only the first possibility appears, the second happens e.g. in the (DN) 
case for the C even breathers). 

The perfect agreement between the TCSA data and the predictions coming from (|6.3| ) 
for B 1 , B 2 and B 3 with /„ > in both the (DN) and the (NN) cases - as summarized in 
Figs.( |6T3| - |6.4D - gives convincing evidence for the correctness of the reflection factors given in 

!■ □ 



6 Among the breather lines these are of the lowest in energy above the ground state. 



7 On Figs. (6^-6^) the dimensionless energy levels above the ground state are plotted against I. On all plots 
the continuous lines are the interpolated TCSA data and the various symbols mark the data corresponding 
to the various BY lines. Some of the higher TCSA lines appear to have been broken, the apparent turning 
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Using the parametric form of L(9) given in ( |6.3|) together with the limits 
hmR^(9)R^(9) = ±1, lim | log R { ^ (9) R^(9)\ < oc 

a— >0 O^oo 

it is straightforward to show that for I n > the range of L(9) runs from zero to infinity, as 9 
decreases from infinity to zero (i.e. as we move from the UV to the IR). However, as it was 
explained in section 4, in the case of I n — there is a maximal value of L beyond which 
we can not go keeping 9 real. We can see this phenomenon on the I n = lines on Figs. (|6.3|- 
IfTil) . Interestingly, there are TCSA lines on Figs.( |67^ - |6^ ) which, below a maximal L, can be 
described well by breather (BY) lines with negative quantum numbers like the B 2 line with 
I2 = —1/2 or the L> 3 line with J 3 = —1. The essential difference between the L{9) functions 
with I n = and the ones with I n < is that while the former ones reach their maximal value 
at 9 = 0, the latter ones do this at some positive 9. Therefore, while in the former case the 
BY lines can be continued by going to purely imaginary rapidity, we can not do this in the 
latter cases. As illustrated on Figs.( |6.3| - |6.5|) this continuation describes the TCSA data very 
well as u = —i9 moves from zero towards the first singularity in the reflection factors, if this 
first singularity corresponds to a bound state pole in bootstrap: see the B 1 lines ending in the 
ei bound state in both the (DN) and the (NN) cases. 

These Figures also show that the idea put forward in section 4 to describe (at least for 
large enough L-s) the TCSA lines corresponding to BBS by appropriate one particle (here 
breather) BY lines with purely imaginary rapidities (independently whether they are obtained 
from continuation through 9 = or not) works indeed. At this point the comparison of the 
DN and NN spectra is very instructive: while in the DN cases there are only single lines that 
descend to t\ or e 2 from above for large L-s, in the NN cases one can see both an increasing 
and a decreasing line tending to t\ and e 2 . In the latter case these nearly degenerate TCSA 
states may be interpreted as the symmetric and antisymmetric wave functions in which one of 
the identical (and hence indistinguishable) boundaries is in an excited state and the other one 
is in the ground state. For large enough L-s the volume dependence of the nearly degenerate 
energies is well described by BY lines with purely imaginary rapidities just above or just below 
the values corresponding to e\ or e 2 . 

On the last two plots in Fig. (|6.5| ) we also exhibit the appearance of the |3) state in an 
expected way. For p = 25/137 (when the state |3) is absent) the BY line marked by the o-s 
is a solitonic line with N = 1 extending to infinity, while for p = 50/391 (when |3) is already 
present) it is a solitonic line with N = 0, ending at a maximal I. In this latter case the TCSA 
data in the continuation of this BY line are correctly described by the x-s, which are the data 
from the BY line of the third breather with imaginary rapidity just below 

6.2.2 Soliton lines 

As mentioned in the introduction to this section in the description of solitonic one particle 
states we have to face the problem of non diagonal reflections. In the C symmetric case (i.e. 
when s «-> s is a symmetry) the one particle soliton anti soliton transfer matrix is 2 x 2, and 
the Bethe Yang equations using the eigenvalues of this matrix take the form 

1 = ^mls^o { Pl{6) ± Q L (e)){P R (0) ± Q R (6)) , (6.4) 

points are in fact level crossings with the other line not shown. This happens because our numerical routine, 
instead of giving the eigenvalues of the Hamiltonian in increasing order at each value of /, fixes their order at 
a particular small I and follows them - keeping their order - according to some criteria as I is changing to 
higher values. 
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;;ure 6.3: TCSA data and soliton/breather Bethe Yang lines in the C even sectors with 
= 25/137. 
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where, in the C even/odd sectors the upper/lower signs apply, and Pl,r(9) (Ql,r{9)) denote 
the soliton soliton (soliton anti soliton) reflection factors on the ground state left and right 
boundaries. In the (NN) case these equations simplify to 

1 = e i2lsinhe {P{6) ± Q{9)) 2 , I = ML, (6.5) 

while in the (DN) case - when on the Dirichlet wall Q vanishes - to 

1 = e i2lsinhe (P(9) ± Q{9))P D {6). (6.6) 

(Here Pd is the £ = case of the soliton reflection factor in BSG with Dirichlet boundary 
conditions: Pd{9) = P { - + \9)\^=q of ref. (4|]). Using the infinite product form of P, Q and Pd 
it is straightforward to verify, that the expressions multiplying e l2 ' smhe in Eq. (|6.5| - |6.6|) are 
indeed pure phases - at least for < 9 < oo - thus these equations make sense. 

The infinite product representation is not very useful when we compare BY lines coming 
from Eq. ( |6.5| - |6^ ) to the TCSA data. Therefore, using the well known integral representation 
for the logarithm of the Gamma function, we recast 1(9) such that the parametric form of the 
solitonic BY lines looks like (with E(l)/M denoting the energy above the ground state): 



j/0x = Nn + ai\ogr± + I(p,9) 



(E(l)/M,l) = (cosh #,/(#)), (6.7) 

where 

1(9) = Nn + m 

sinh 9 
with 

oo 

J(p, 6)=f^ sin(^) ( 2smh(3y/2)sinh([l-p]y/2) + asinh^) - sinn^ x 
J y 7r V sinh(py/2) sinh(2?/) cosh(y) sinhyo?/) / 

o 

and 

r + = r_ = (6.10) 



in the C even/C odd sectors. Here the parameter a is one, a = 1, in the (NN) case and is 
one half, a = 1/2, in the (DN) case; the quantum number N, characterizing the BY line, is a 
non negative integer in all sectors and all cases and -as before - 9 is both the parameter of the 
line and the rapidity of the solitonic particle. (Although this representation of the solitonic 
BY lines works for real rapidities only, one can show that all the "boundary dependent" poles 
at 9 = iu come from logr-t, in fact from the set of poles, 9 = z(| — kpn), r + has poles for 
k = 2N while r_ for k = 2N + 1). 



On Figs. (|6.3[ - |6.4| ) we demonstrate the excellent agreement between the TCSA data and 
the BY lines given by ( |6.7| - |6.10| ) with > in both the C even and the C odd sectors of 
the (DN) and (NN) cases. This agreement gives a strong evidence for the correctness of the 
solitonic reflection factors given in Q. 

Note that in some cases we have TCSA lines that below a maximal L can be described well 
by solitonic BY lines with zero quantum number, = 0. The continuation of these BY lines 
to imaginary rapidities - that worked for the breather lines - is hampered by the fact, that in 
all the cases investigated, the pole in (P±Q) 2 or (P±Q)P D nearest to the origin turned out to 
be a Coleman-Thun pole. Therefore the solitonic BY lines continued to imaginary rapidities 
deviate from the TCSA data. 
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gure 6.4: TCSA data and soliton/breather Bethe Yang lines in the C odd sectors with 
= 50/391. 



24 



6.3 New boundary bound states and reflections on excited walls 

The most interesting application of TCSA is to find evidence for the existence of the new 
BBS-s predicted by the bootstrap and check the correctness of the 'excited' reflection factors, 
i.e. the ones describing the scattering of breathers/solitons on the wall in an excited state. In 
this inquiry it is very helpful to analyze the DN spectra: if, e.g. we find a TCSA line (state) 
with asymptotic (large I) energy e\ + e 2 in the DN spectrum, then we can be sure that it 
indeed corresponds to the |1,2) state of the Neumann boundary. (Note that in the NN case 
there are more states having this asymptotics, as the energy of the configurations, when one 
of the Neumann boundaries is in the |1) state and the other is in the |2) state, also tends to 
ei + e 2 ). 

The state |1,2) can be generated by the third breather, furthermore, the corresponding 
pole at Mi is the one nearest to the origin in the third breather's reflection factor. Therefore 
it is expected that the state |1,2) becomes visible as the imaginary continuation of a B 3 BY 
line with J 3 = 0. We collected the evidence for the existence of |1,2) on fig. (|6.5| ). On all 
three plots the continuous lines are the interpolated TCSA data, the +-s and the empty boxes 
denote the BY line of B 1 and B 3 with ii )3 = respectively, while the full triangles and the *-s 
denote the continuations of these lines to purely imaginary rapidities with u ranging from zero 
to u\. Note, that for p = 2/7, when, according to bootstrap considerations, there is no |1,2) 
state and the imaginary continuation of the B 3 BY line should run into a Coleman - Thun 
pole, the BY line indeed departs from the TCSA data, (which, in this case are well described 
by a B 2 B l BY line where the second breather's rapidity is imaginary and the first breather's 
quantum number is I\ = 1, the data corresponding to this line are given by the m-s). On the 
other hand, when p is such that |1, 2) should exist, we see that the imaginary continuation of 
the third breather's BY line describes the TCSA data just as well as that of the first breather's 
line, which runs into e\. 

Further evidence for the existence of a new state predicted by the bootstrap is given on 
Fig. (|6.6|) , where we exhibit a single TCSA line descending to t\ + e 3 . In the light of what we 
said above this indicates the existence of the boundary bound state |1,3). (The horizontal 
lines corresponding to m 4 = M 4 /M and the dimensionless t\ + e 3 are drawn to guide the 
reader's eye to distinguish this line from the real fourth breather's BY lines lying above m 4 .) 
This state corresponds to the pole at u 2 in the fourth breather's reflection factor on the ground 
state boundary. Therefore, above a certain minimal I, the TCSA data in this line can certainly 
be described as the BY line of B^' with imaginary rapidity below U2, but, since this reflection 
factor also has a pole at U\, can not be obtained as the analitycal continuation of a BY line 
through = 0. 

The bootstrap procedure gave not only the spectrum of the BBS but also the reflection 
factors describing the scatterings of the various breathers/solitons on excited boundaries. On 
Fig. (|6.6|) we also confirm the correctness of these new reflection factors. 

There are two ways one can describe the scattering of a breather, (B n say), on an excited 
boundary \k). The first is to use a BY line ("excited" BY line) following from Eq. (|4.4j) 

(E(L),L) = (e k + M n coshe,(27rI n + i\ogRP(6) +ihgR^\e))/(2M n smhe)) , (6.11) 

where or R^ now corresponds to the appropriate Rj?) , and the energy above the ground 
state is given by e k + M n cosh9. This description, however, ignores the fact, that, in finite 
volume L, even in the absence of the moving B n , the energy of the state \k) only asymptotically 
coincides with e k . To remedy this we may try to use a two particle BY line, where ground 
state reflection factors appear, but one of the particles moves with imaginary rapidity, (which 
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p = 50/391, the |1,2) state exists 
Figure 6.5: Evidence for the existence of the |1,2) state from three C = — 1 DN spectra. 
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Figure 6.6: TCSA data, the state |1,3), and "excited" BY lines in the C even sector of the 
DN model with p = 50/391. 

is in the vicinity of a pole corresponding to \k) in the appropriate reflection factor). It is 
expected that in this way we describe more accurately the volume dependence of the energy 
of the state combined from the moving breather and \k). 

Since the TCSA data are usually the better the lower the level is, we may hope to confirm 
the new reflection factors mainly for light particles reflecting on lowly excited boundaries. In 
the p range we investigated the first breather is the lightest particle, therefore on Fig. (|6.6|) 
we compare B lJ s BY line on the first excited boundary |1), and the TCSA data. The points 
marked by * correspond to 1\ — 1, 2 in ( |6.11| ). The data, corresponding to the alternative 



description, using two particle B 1 B 1 lines with one imaginary rapidity and one real one with 
quantum numbers I\ — 1 . . . 3 are marked by the full triangles; they give a somewhat better 
description of the TCSA lines. The four lines marked by the empty boxes correspond to two 
particle B 1 B 1 lines with real rapidity and quantum numbers (1,0), (—1, 1), (1,2) and (—2,2) 
respectively. The lines marked by o and + are soliton with N = 1 and B 2 on ground state 
boundary with = |, and = |, respectively; their function is to help the reader to 
identify the already known breather/soliton ground state BY lines. Finally, for comparison, 
we also included in this plot two lines - marked by «-s; they are solitonic "excited" BY lines 
with N = 1,2. 

Summarizing, we can say, that for the first breather, the satisfactory agreement we find 
between the TCSA data and the various ways of describing its scattering on the first excited 
boundary gives support for the correctness of the new reflection factors. We expect that the 
relatively poor agreement between the (considerably higher) TCSA data and the "excited" 
solitonic BY lines could be improved by replacing the latter by two particle soliton- .B 1 lines 
where the breather is moving with an imaginary rapidity. On Fig. ( |6.6j ) we also show examples 
of real two particle BY lines correctly describing TCSA data as well as instances, when a TCSA 
line is described by several BY lines. In the C odd sector the lowest lying '"excited" BY lines 
describe the scattering of B 2 on |1); we repeated also in this case the analysis described above 
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and obtained the same qualitative conclusions. 



7 Conclusions 

In this paper we analyzed the sine-Gordon model with Neumann boundary condition (SGN). In 
particular we established the spectrum of boundary states, developed a framework to describe 
finite size effects in boundary theories and used this framework together with TCSA to confirm 
the boundary states and the reflection factors. 

Of the boundary states we showed that they can be labelled by an increasing sequence 
of positive integers satisfying rii < [|1 , and denoted as . . . n k ). For the states |n) 
we showed that their energies coincide with that of the WKB quantized classical boundary 
breather. We realized the existence of the states \n, k) by studying on the one hand the poles 
of the breather's reflection factors on non excited boundaries, and on the other, the poles of 
the soliton reflection factor P\ n )(u) on the excited state \n). The novel feature of this study was 
that we computed also the residues of the Coleman-Thun diagrams and the reflection factors. 
This way we discovered that Coleman-Thun diagrams and creation of boundary bound states 
may coexist: this happens if the residue of the diagram is not sufficient to explain the residue 
of the reflection factor. Finally we proved that the new boundary states labelled by more 
than one integers can be generated successively by appropriate soliton reflections. We also 
determined all solitonic and breather reflection factors on these new states and explained their 
poles. 

It is interesting to compare the labeling of boundary states here in SGN to that of in 
sine-Gordon model with Dirichlet boundary conditions |§]. In the latter case the labeling is 
more complicated: it consists of two sequences of positive integers. The explanation is that in 
the Dirichlet case, in contrast to the Neumann one, the solitons and the anti solitons generate 
different boundary states. 

Moving from the infinite half line to a finite line segment, we gave the general form of 
Bethe-Yang equations for particles in finite volume with boundaries and used them to discuss 
the appearance of boundary bound states in this finite size setting. We showed in particular, 
that these states can be described as solutions of the BY equations with purely imaginary 
rapidity. 

Finally we studied the finite volume spectrum of SGN by the truncated conformal space 
approach. We established, that the ground state energy density of this model coincides with 
that of the bulk SG and could also relate the ground state Neumann boundary energy to that 
of the Dirichlet one with 4>q = 0. The perfect agreement we found between certain lines in 
the TCSA spectrum and the one particle soliton/breather Bethe-Yang lines gives evidence 
that the reflection factors entering into the BY lines are indeed correct. We showed that the 
boundary bound states indeed appear in the TCSA spectra as predicted by the Bethe-Yang 
equations with purely imaginary rapidity. Finally in the TCSA data we found evidence for 
the existence of the new multi labelled states and also for the excited state reflections. 
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Figure A.l: Schematic draw of the top line 



A The proof of the existence of boundary bound states 

In this appendix we prove that the boundary bound state \ni, n 2 , . . . , ri[, k) can be created 
by reflecting a soliton with rapidity v k on the boundary state |n 1? n 2 , ■ ■ ■ , ni). In doing so we 
show that no on-shell diagram exist for this process. Since any on-shell diagram consists of 
virtual processes we omit the attribute virtual from now on. One possible diagram is shown 
on Fig. O. 



The leftmost process in which a soliton arrives with rapidity and virtually fuses with 
a breather into a soliton is called the final process while the breather is the final breather. 
Now if we follow the breather line towards the wall and turn in each 3-point vertex to the 
left we obtain a line which we call the top line. (In case of a 4-point vertex, that describes a 
possible two particle scattering instead of fusion, we follow the line with the same rapidity). 
The top line reaches the wall in the initial breatherf]. Of course we have analogous diagrams 
by replacing some of the breather lines with soliton lines. In this case the final/initial soliton 
terminology is used. Clearly in the final process the particle created must travel towards the 
wall i.e. it must have positive rapidity. 

The outline of the proof goes as follows: First of all we introduce the notion of "proper" 
particles and show that particles with rapidity suitable for being final ones are "proper". Since 
any initial particles are also "proper" we show that in the diagram above only "proper" particles 
exist. Moreover, we determine all the bulk processes that take place on the top line and call 
them "proper" processes. Having introduced a charge like quantity for "proper" particles we 
show that the maximal charge of the initial particles is n\. Since the charge never increases in 
"proper" processes and the charge of the final particle is at least k which is larger than ni we 
conclude that no on-shell diagram exists. 

Let us go into the details. We call a breather "proper" if it has rapidity — u n , n & N and 
a soliton is "proper" if its rapidity is —v n . The possible final processes are shown by diagram 
(a-b) on Fig. (|A.2| ). Clearly both the final breather and the final solitons are "proper". Initial 
particles are produced by the decay of boundary bound states, the various virtual processes are 
described on the three diagrams of Fig.([Aj|). Clearly all the initial particles are also "proper". 

8 We call it in this way, because in our convention, when time flows from top to bottom, the emission of 
this breather by the wall precedes all the other processes on the top line. 
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Figure A. 2: Final processes 
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(b) Breather initial state 



(c) Soliton initial state 



Figure A. 3: Initial processes 



For "proper" particles we introduce the notion of charge. For B n with rapidity u 



charged is defined to be Q 



n+m 



-u m the 

while for a soliton with rapidity u = — v k it is simply Q = k. 
The charges of the initial particles are also displayed on Fig. (|A.3| ) . The maximum of the initial 
charge is ni, obviously. The final charges are also displayed on Fig. ([A.2|) . Since the particle, 
produced in the final process, has to travel towards the wall, (i.e. v n _ k < v k for diagram (b)), 
we conclude that the final charge must be larger then k. Now we analyze how the charges 
change on the top line. 

We are looking for bulk processes which can take place on the top line i.e. which are in 
between the final and initial "proper" particles. There are two types of bulk processes, one 
with only breathers and the other one with two solitons and one breather, already depicted 
on Fig.( [3.1|) . In the present case however, any of the particles can arrive from the boundary 
so we have six different cases. First we analyze how the rapidities are changing. Consider the 
diagrams on Fig.( |A.4| - |A.5|) . 

We start with the case when the final particle is a breather. Let us denote, for a moment, 
by N x the number of x type process on the top line, where x runs through the six diagrams 
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Figure A. 4 



on Fig.( [A.4|) and Fig. (|A.5| ). From the conservation of particle types it follows that in case 



of a breather initial state ./Vb rea ther- soliton — A^nton-breather, while in case of a solitonic initial 
state iVbreather-soiiton + 1 = ^soiiton-breather and all other N x s are unconstrained. We claim, 



however, that only diagrams (a-b) on Fig. (|A.4|) are allowed. In showing this we parameterize 



the rapidities as fT| + u n and note that for a "proper" breather K = 0, while for a "proper" 
soliton K — — 1. Now it is easy to see that on diagrams (a-b) on Fig.( |A.4| ) and on diagram (c) 
on Fig.( [A.5| ) the parameter K does not change, on diagram (c) on Fig. (|A.4| ) K — > K + 2, on 
diagram (a) on Fig.( |A.5| ) K — > K — 1, finally on diagram (b) on Fig.( |A.5| ) K — * K + 3. Since 
the initial and final particles are "proper" the various i^-factors must sum up accordingly 
on the top line, showing that only diagrams (a-b) on Fig. (|A.4|) are allowed. We call these 
processes "proper". 

We also indicated on the diagrams how the charge changes. If we did not have diagram 
(b) on Fig. (|A.4j ) we would be ready, since in the process on diagram (a) the charge never 
increases, consequently the final charge can not be larger than the initial one. In the general 
case, having also the breather fusion type process we argue as follows: find the breather fusion 
vertex on the top line which is farthest away from the wall. The bottom breather coming from 
the direction of the wall in this process is also a "proper" breather. Following its line back 
towards the wall always turning to the left we either reach the top line or the wall. Thus this 
new line has one or two common vertices with the top line. Now we change from the original 
top line to a new "top" line by replacing the segment of the the old one between the common 
vertices or between the common vertex and the wall with the line obtained from the bottom 
breather. Since the new line segment shares the properties of the original top line, namely it 
starts and ends with a "proper" breather and we have turned at each vertex to the left, we 
can argue for it in the same way as before. Applying this procedure step by step we arrive 
at a line which contains only breather decay type processes so the previous argument applies. 
Since in every step we eliminate at least one edge of the diagram the procedure terminates and 
we arrive at the final breather concluding that its charge can not be larger than the largest 
initial charge i.e. rii. 

For the solitonic final state the K parameter is —1. Thus either we have purely solitonic 
processes on the top line, in which case the charge decreases with each step or we start with 
some breather processes allowed in the above mentioned sense and later by the breather- 
soliton diagram on Fig.( [A.5|) we change back from the breather to soliton and continue with 
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(a) Breather-soliton 



(b) Soliton-breather 



(c) Soliton-soliton 



Figure A. 5 

the purely solitonic processes. Now in the second case the previous proof for the breathers 
can be used to show that the charge does not increase in the breather part of the top line. In 
the change from a breather to soliton the charge always decreases, since for B n with charge k 
the inequality n < 2k is always satisfied, so the conclusion is the same i.e. the final charge is 
always less than the maximal initial charge rii. 

Summarizing we showed that no on-shell diagram exists, since for the existence the final 
charge must be at least k, but the initial charge is at most n\ < k and the charge does not 
increase in the allowed processes. 

These considerations apply to the direct channel process only. The proof for the crossed 
channel is quite similar so we omit it here. The difference is just that we have to introduce a 
new charge for B n which is q = if the rapidity is u = u m . Having analyzed the analogous 
diagrams step by step the conclusion is the same. 

B Pole structure 

In this Appendix we start by reviewing the soliton and breather reflection factors on the state 
|ni,7i2, • • • ,nk) and analyze their pole structure. Proceeding inductively in the order of the 
pole we show nested on shell diagrams which have poles of the same order as the reflection 
amplitudes. We do not compute the residue of these diagrams, however, since on the one hand 
the calculation is very difficult to perform, and on the other, we know instances where the 
result does not coincide with the one coming from the reflection factor. The analysis is carried 
out in two steps. First, the solitonic poles are explained, then we turn to the explanation of 
the breather ones. 

B.l Solitonic pole structure 

The solitonic reflection factors are 

P\ ni ,n 2 ,...,n k ){u) = a ni (u)a n2 {u) . . .a nk (u)P{u) ; P^Q. 

The function a n (u) = {2n — 1 + \}{2n — 3 + A} ... {1 + A} has the pole structure shown on 
Fig. flB.lj ), where the number of dots represents the order of the pole. 
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Figure B.l: Poles of a n {u) 



In the general reflection amplitude we have to sum up all the orders of the poles of the 
various a ni (u) factors and also to take into account the poles of the prefactor P(u). As a 
result, we have the following table 
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Now the main observation, similar to what was made in |4]], is that the poles above are 
already present in the reflection amplitude P\ nk )(u), only the orders are different, but this can 
be compensated by nesting appropriate sub-diagrams as we will see. Let us consider the poles 
at u N step by step. 



• As we proved in the previous Appendix the simple pole for N > n k is responsible for 
the creation of the state \ni, n 2 , . . . , n k , N). 

• Lets consider diagram (a) on Fig. QB.3; ). The state \ni,n 2 , . . . ,n k _i,n k ) decays to 

\rii, n 2 , . . . , n k _i) by emitting an anti soliton with rapidity u nh . The soliton and the anti 
soliton fuse into B N+nk with rapidity u nk -N- In the N = n k case the breather travels 
parallel with the wall giving a diagram with a second order pole. 

• In the N < n k case the reflection factor of B N+7lk has a pole at u — u nfe _7v, for which 
a sub-diagram, that is described in detail in the next subsection, has to be embedded. 
We just present here in advance, that when N decreasingly reaches any of the n^-s then 
the order of the breather pole is increased by one and when N = rii — 1 it increases 
by one more. Furthermore, it does not change until N reaches rij_i, where the same 
story happens. Collecting all the order of these poles and also the one more extra, which 
comes from the ground state reflection factor we arrive at the right result. 

We emphasize that from the shifting argument described in section (3.3) we know that in spite 
of the existence of the on-shell diagram shown above the pole at in the reflection amplitude 
P\ni,n2,-,n k )( u ) a ^ so creates the state |ni, . . . ,rii, N,n i+ i, . . . ,n k ), if this state exists, that is if 
rii < N < n i+ i for some i. 



B.2 Breather pole structure 

The breather B n on the state \rii,n 2 , . . . ,n k ) has the following reflection factor 
where 

un i \ r ,o ii r io qi r io 1 on / {n - 2m + 1} if - 2m + 1 > 
b m {u) = {n+2m-l}{n+2m-3\ . . .{n+2m-l-2l\ . . . < } { ,i^n 

m L J [ |2m — n + 1} if 2m — n + 1 > 
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g n+x-np 



g np-x sjip+x 



In ,,..,n ...,n . > 

1 ' p k 



In ,,..,n+x,..,n .> 
1 ' k 



ln r ..,n r .,n k > 



In ,,..,n ...,n , > 




(a) Breather decay 



(b) Soliton decay 



Figure B.2 



Lets parameterize the poles of U^iu) as u n+2x - Considering the two cases separately we obtain 
the location of the poles as 



o o 
o o o 



o 
o 



o 
o 



-m 



in the case when n > 2m while as 



o o 
o o o 



o 
o 



in the other case. Summing up all the poles coming from the various b -factors we get 
the following result. If n q < x < n q+ i and n p < n + x < n p+ i hold in the x > case or 
n q < M < n q+i and n p < n — \x\ < n p+ i is satisfied for x < 0, then the order of the pole at 
u n+2x is 2(p — q) + e (where e = 1 if the bound for p is saturated and e = —1 if the bound 
for q is saturated and zero otherwise). For each pole there is the analogous pole at tc — u n+ 2 Xl 
which is in the physical strip if u n+ 2 X is not there. 

Now we explain the pole of the reflection amplitude R^ n2 nk )( u ) a ^ u n+2x- For x = n k 
the state |ni, ■ ■ ■ , nk + n) is created or if it is not in the spectrum we have an analogous 
diagram to diagram (d) on Fig. (|3.2j) . For x < n k consider diagram (a) on FiR. (|B.2|) . 



Here B n decays into B n+X ~ np and B np ~ x in the bulk. First B n+x ~ Up reaches the wall 
with rapidity u n+x+rip exciting the state \n\, . . . ,n p , . . . , n^) to \n\, . . . , n + x, . . . , n^). If the 
breather traveled backwards in time (n + x + n p > [A]) or the excited state did not exists 
(n + x > [!]) then we would have diagram (b), where the state \ni, . . . ,n p _i,n p ,n p+ i, . . . ,n k ) 



would decay to \ni, 



n 



n k ) by the emission of a soliton with rapidity — v np , 



which then absorbs B n+X ~ np and reaches the wall at rapidity — v n+x . In both cases B rip ~ x 
reflects on the resulting state with rapidity U( np - X ) + 2 X - In order to compute the order of its 
pole we observe that n q < x < n q+ i as before, but now n p _i < {n p — x)+x. From this it follows 
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(a) Soliton decay (b) Soliton decay crossed (c) Breather decay crossed 

Figure B.3 

that the order of the pole is 2 (p — 1 — q) + e, which is less than the one investigated originally by- 
two. Since the diagram itself gives a second order pole (we have used the Coleman-Thun type 
cancellation in the (b) case) we are ready. Resolving inductively the pole of this reflection by 
the analogous sub-diagram we arrive at cascades of diagrams, where in the last sub-diagram 
the breather reflects with a simple pole or without poles. 

For the poles of type n — u n+2 x we have an analogous consideration: If we parameterize 
the rapidity as n — u n+2 x and analyze the relation between x, n + x and n.j we get the same 
formula for the order of the poles as we had before. The explanation is given in terms of the 
crossed diagrams. Diagram (b) on Fig.( |B.3| ) explains the x = n k case, while for x < n k we have 
diagram (c) on the same Figure. The argument is completely analogous to the non-crossed 
version so we omit it here. 

Note, that these breather diagrams are generic, which means that they exist for any x. It 
may happen, however, that for some particular x we have some other diagram with the right 
order of pole. Notice also that in spite of the existence of the diagrams above we might have 
boundary bound state creation, described in detail in section (3.3), at the same time. 

C Boundary c = 1 theories 

In this Appendix we warm up by considering the compactified free boson with periodic bound- 
ary conditions on a circle of circumference L. Having mapped the system onto the plane we 
identify those conformally invariant boundary conditions which also preserve the underlying 
affine U\ X U\ symmetry. They originate from cutting the original circle and applying Neu- 
mann or Dirichlet boundary conditions on both ends of the strip. We also determine the 
theory which corresponds to imposing Neumann boundary condition on one end and Dirich- 
let boundary condition on the other end of the strip. In all cases the data we need for the 
TCSA, such as the spectrum of bulk and boundary primary operators, their normalizations, 
conformal weights, bulk and boundary operator product expansions are summarized. Though 
most of these data are available in the literature we collect them here to make the paper self 



35 



contained. 



C.l Compactified free boson with periodic boundary conditions 

Consider a free boson with compactification radius r confined into a box of size L subject 
to periodic boundary condition $(L, £) = $(0, t) + 2nrn, where n G Z is called the winding 
number. The action which governs the dynamics is 



S7T 



S = — I dt I dxd^d^ 



-oo JO 

Canonical quantization results in the following expression for the field operator t): 



t) = $ + — ( U t + —x j + i 2^ - (a n e 1 ^"^ + a n e 



where M is the winding operator with eigenvalues m G Z, and the nonzero commutators are 

[$o,n ]=2 ; [a n , a m ] = n5 n+m , [a n , d m ] = n5 n+m . (C.l) 

We can use (^-function regularization for computing the Casimir energy and find the Hamilto- 
nian as 

Sttv^ _ 1. 2tt /, frM\ 2 \ 

H = y^( fl - A + ' I - A " 12' + T I 0+ V~2~/ /' 

Having mapped the system onto the plane via the (x, it) = £ — > 2 = e*^ conformal transfor- 
mation, chiral factorization originates from the split §(z,z) = (p(z) + (p(z), where 

4>{z) = 0o — iao hi z + i y a n ; (j)(z) = (f) — id In z + i } a n . 

^— ' n z — ' n 

It is useful to introduce the dual of $(z, z) as $(2;, z) = <p{z) — <p{z). The U\ x Ux symmetry 
of the model is generated by the chiral currents, J(z) = id z (f>(z) , J{z) = id 2 (j)(z), and the 
primary fields of the symmetry algebra are the vertex operators 

V {n>m) (z, z) = : e ^^-z)^H^z) , = , e i q m+m*) : (c.2) 

We use the parameterization (n, m) and (q, q) in parallel, the connection between them being 
<? + <? — an d q — q = rnr. The vertex operators have conformal weights /i( n>m ) = \ , h( n ,m) = 
^- with respect to the conformal energy momentum tensor T(z) = | : J(z)J(z) : and = 
I : J(z)J(z) :. The Hilbert space is built up by the successive application of the modes of the 
chiral currents on the highest weight state created by the vertex operators as 

\n,m) = V M (0,0)| 0) 

and ilo|n, m) = n\n, m) , M\n, m) = m\n,m). The operator product expansion of the vertex 
operators starts as 

W^Vfc^KtD) = {z-w)«\z-w)rfv M j Mf) (w,w) + ... (C.3) 



We have two maximal sets of the allowed (n, m) pairs QT9|| and here we concentrate on the 
bosonic one which has n G Z, m G Z. 
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C.2 Conformal boundary conditions 

Now we would like to impose such boundary conditions in the theory that preserve not only 
the conformal but also the Kac-Moody symmetry of the model. The easiest way to formulate 



this is to restrict the theory on the upper half plane and demandpO 

L(z) = L(z) , J(z) = fij(z) ; 



on the real axis, where Q is any automorphism of the algebra (|C.1|) . We have two possible 



choices Q = 1 ov Q = —1, plf , they correspond to Neumann and Dirichlet boundary conditions 



imposed on the field <3>, respectively. Note that for the dual field the roles are interchanged, 
that is Q = 1 corresponds to Dirichlet while Q = — 1 for Neumann boundary condition. 

C.2.1 Neumann boundary condition 

The Q = 1 boundary condition in terms of the modes reads as a n = a n , that is 

$(z, z) = $ - iooOn z + In z) + i V — (z~ n + z~ n ) . 

z — ' n 

We can map this field back to the strip by the inverse of the transformation £ — > z = e*^ 

t ) = $ + ^n + V ^L( e *f + e - J r «(-+*)) (C.4) 

From the x -dependence we can read off the Neumann boundary conditions 

d^{0,t)=d^{L,t) = ; Vt 



The Hamiltonian of this system in the strip can be obtained directly from (|C.4|) by C-fimction 
regularization 

= — + - JJna_ n a„ - — . (C.5) 

Turning back to the plane we realize that the vertex operators V( n ,o){z,z) do not reproduce 
the most singular part of the OPE ( |C.3|) any more as we would expect. Changing however 
their normalization function as 

V (n , 0) (z, z) = \z- zf ■ e iq ^ {z ~ z) ■ ! Q = ~ , (C6) 

2 

we preserve their conformal weights, (which is h — %- for T(z)) and restore the leading part 



of the OPE (|C.3|) . The boundary fields in a boundary conformal field theory live only on the 
boundary and are in a one-to-one correspondence with the vectors of the Hilbert space. In 
our case they are the boundary vertex operators 

M'„(.r) ; q = 

r 

(which are primary fields of weight h n = 2q 2 ), and their descendants. The dual field vanishes 
at the boundary so in the Hilbert space we do not have winding as we would expect from the 
Neumann boundary condition. 
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C.2.2 Dirichlet boundary condition 

The Q = — 1 boundary condition forces $(z, z) to be of the form 



${z, z) = %- ia (\nz - In z) + % V — (z~ n - z~ n ) 



The boundary condition can be determined on the strip by mapping the system via the inverse 
of £ — > 2; = e*i^ to the strip. We have 

L n ' 

where M = M + and M G Z is the winding number. This corresponds to the 

$(0,t) = $ ; $CM)=$ L Vt 
Dirichlet boundary condition. The Hamiltonian of the system above is 

2tt / rM\ tt /v-^ 1 



By changing the normalization function of the vertex operators 



V m (z, z) = \z- z\- q2 ■■ e iq *^ : ; q = £ , (C.8) 

we can ensure the most singular part of the bulk OPE to hold. Now the field <3> is constant at 
the boundary so we do not have nonzero momentum Fock modules. The dual field however 
can create the sectors corresponding to the different winding numbers. The Hilbert space by 
this token consists of the Fock modules corresponding to the different winding numbers. 

C.2.3 Mixed boundary condition 

It is also possible to demand mixed boundary conditions. On the plane one possible choice 
is Dirichlet boundary condition for ^sm(z) = , ^Ste(z) > and Neumann for ^sm(z) = 
, $le(z) < 0. A field satisfying this can be given as 



^(z,z) = % + t V ^( z - n -z- n ). 

n 



Since on the strip 



n 



the boundary condition there reads as 

$(0,t)=$„ ; d x $(L,t) = Vf, 



38 



that is Dirichlet boundary condition at x = and Neumann boundary condition at x = L. 
Using ^-function regularization the Hamiltonian turns out to be 



I Yl na ~n a n + J • (C.9) 



The model looks like a twisted boson with chiral current J(z) = Yl a n z n 1 • The vertex 
operators acquire a factor 

V q (z,z)= ( .^V' 1 ; g = -, (C.IO) 

in order to ensure the most singular part of the bulk OPE to hold. Now the field $ has a 
non vanishing limit only for 3m(z) = , 9?e(z) < while the dual field survives at ^sm(z) = 
, 9fte(z) > 0. Thus the only boundary operator is the current J(x) and consequently the 
Hilbert space is generated by acting with its modes a_ n , n £ Z + — |, on the vacuum vector. 
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